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Abstract We demonstrate that the fractional Laplacian (FL) is the principal characteristic operator of 
harmonic systems with self-similar interparticle interactions. We show that the FL represents the “ fractional 
continuum limit’ of a discrete “self-similar Laplacian” which is obtained by Hamilton’s variational principle 
from a discrete spring model. We deduce from generalized self-similar elastic potentials regular representations 
for the FL which involve convolutions of symmetric finite difference operators of even orders extending the 
standard representation of the FL. Further we deduce a regularized representation for the FL —(—A) 2 holding 
for a £ IR > 0. We give an explicit proof that the regularized representation of the FL gives for integer powers 
j £ No a distributional representation of the standard Laplacian operator A including the trivial unity 
operator for a —> 0. We demonstrate that self-similar harmonic systems are all governed in a distributional 
sense by this regularized representation of the FL which therefore can be conceived as characteristic footprint 
of self-similarity. 
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1 Introduction 


Fractional calculus has a long history, but it took until the last decade for a new interest to emerge to employ 
fractional operators and the so called fractional Laplacian (FL) (often also referred to as Riesz fractional 
derivative) — (— A) 2 where E IR indicates a fractional, in general non-integer, power of the standard Lapla¬ 
cian A where the sign is chosen such that the FL is for a > 0 a negative (semi-) definite self-adjoint operator, 
formally defined by its Fourier representation — k a < 0. The positive exponent a > 0 is a requirement which 
comes into play when translational invariance (uniform translations corresponding to the k = 0 mode) give a 
a vanishing contribution. The Green’s function kernel of its inverse operator — (—A) - 2 5 the so-called Riesz 
potential, is also well known in the literature and was first introduced by the Hungarian mathematician Marcel 
Riesz m- 

This newly emerging interest in fractional calculus is due to the fact that it has been recognized that 
critical and complex phenomena with self-similarity properties can be often described by fractional equations 
of motion; especially phenomena of anomalous diffusion (Levy flights) are governed by fractional diffusion 
equations. 

One incentive to “popularize” fractional calculus in a broader physics and engineering community was 
probably the survey article of Metzler and Klafter [23] with a broad overview on applications demonstrating 
the power of the fractional approach to model complex phenomena in various fields of physics and especially the 
interlink between self-similar processes such as for instance phenomena of anomalous diffusion (Levy flights) 
and its fractional continuum description. Prominent examples of Levy flights are phenomena as various as the 
time evolution of the stock market [8] and the chaotic particle transport in turbulent flows |20| . Recently many 
models where developed which employ the FL in various physical contexts, among them the description of 
“complex” dynamic phenomena | I S. fHJ [ I2T [ l22 [ l23[ I25 [ [26l 127] I2S1 [29l Ljl and see also the numerous references 
therein. We also would like to mention the recently emerging great impact of fractional Quantum Mechanics 
coined by Laskin m- 

In a recent article we exhibited the interlink between harmonic self-similar spring systems with elastic 
potentials of Weierstrass-Mandelbrot function type and the FL defined rigorously by the “fractional continuum 
limit” m which will be the starting point of the present paper. Here our aim is to demonstrate that the 
FL is a natural consequence of self-similarity (scale-freeness) of elastically harmonic systems. We deduce a 
regularized representation of the FL which holds for any positive exponent a > 0 extending the standard 
representation of the FL which is limited to exponents 0 < a < 2. 

So far, the spatial representations of the FL are limited due to its hypersingular characteristics to 0 < a < 2 
where a — > 2 approaches asymptotically the standard Laplacian A and a — > 0 the unity operator. For a 
mathematical analysis on some key properties of the FL we refer to textbooks such as (Ij. Generally the 
FL is a nonlocal convolutional and highly singular operator. Its standard spatial representation given in the 
literature exists in its non-regularized form only in the interval 0 < a < e.g. [26]. In many physical contexts 
power law spectra and power law distributions indeed are found within this range [3]. However, examples 
exist which justify to analyze the FL also outside the a-range (0,2): Gopikrishnan et al. found a power law 
distribution of fluctuations of financial market indices which exhibit a scaling exponent a being well outside 
the Levy range 0 < a < 2 [5|. These probablity distributions are described by characteristic functions (Fourier 
transform of the probability distributions) given by exponentials ~ exp —k a r (with the scaled time r > 0) 
with scaling exponents a outside the Levy regime (0, 2). Such stochastic processes have to be described by a 
Fokker-Planck equation which employs a spatial representation for the FL which exists outside 0 < a < 2. To 
render such problems accessible to the analysis, it is necessary to have such a spatial representation for the FL 
with an extended a-range of validity, and it is more than ever desirable to have a spatial representation of the 
FL which covers the entire range of principally admissible exponents a > 0. Such a representation for the FL 
outside the Levy range (0, 2), however, requires regularization since the standard representation of the FL, in 
its common unregularized form, does not exist outside of this interval 0 < a < 2 ( more precisely: for a > 2 
the convolutional integral of the unregularized standard representation of the FL diverges as ~ r 2 ~ a at r —> 0). 
The goal of the present paper is hence to deduce by means of a simple regularization rule a regularized unique 
spatial representation for the FL which holds for any a > 0. 

In the literature as well as in the present paper we define the FL such that it has all “good properties “ of 
the elliptic, standard Laplacian (having Fourier transform —k 2 ). These good properties of the FL are: self- 

1 In the stochastic context of Levy statistics the range 0 < a < 2 is refered to as “Levy interval” [§]. 
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adjointness, spatial isotropy and negative-semidefiniteness (“semi-” since the Fourier transform — k a of the 
FL is vanishing at k = 0 reflecting translational invariance). The latter property requires positive exponents 
a > 00 

In the present paper the FL operator is formally defined as the operator —(—A) 2 having negative eigen¬ 
values — k a = —(k 2 ) 2 < 0 for any a £ IR > 0 (vanishing only for k = 0 when a > 0), leading to all above 
mentioned good properties of a Laplacian. However, this formal definition of the FL does not provide any 
spatial representation of this operator. Our main goal is to deduce spatial representations for the FL from 
simple elastic potentials in a rigorous manner by application of Hamilton’s variational principle. We will show 
in the present paper that the FL is the “natural” operator, and in a sense the most basic operator that is 
generated as a continuum limit from a physical “self-similar” spring model and its generalizations. 

The paper is organized as follows: As point of departure in section [2] we evoke a ID harmonic spring model 
with harmonic elastic potential energy which describes self-similar interparticle interactions which we devel¬ 
oped recently [9]. This discrete model leads to fractal dynamic vibrational characteristics such as a dispersion 
relation of the form of Weierstrass-Mandelbrot fractal functions. Then we analyze a “fractional continuum 
limit” which we introduced earlier m where the self-similar Laplacian of the discrete model takes up to a 
normalization factor the form of the standard representation of the FL which is employed in many references, 
e.g. [HI E3 E3 ESI E>UI3] with restricted exponent 0 < a < 2. We demonstrate that the discrete fractal 
model is the natural discrete counterpart of the fractional model. Both models are linked by the “fractional 
continuum limit” which we introduced recently m- Further in section [A2| we deduce spatial representations 
of the FL which hold for an extended range of exponents a. To this end we introduce generalized discrete 
self-similar harmonic potentials which involve higher-order finite differences in the displacement field of integer 
orders m G HI and analyze the fractional continuum limit: Hamilton’s variational principle then defines in a 
rigorous manner self-similar Laplacians leading in the fractional continuum limit to new representations for 
the FL which exist for exponents 0 < a. < 2m where m E I'll indicates arbitrary integers. These representa¬ 
tions depend on the order m £ W of the finite difference operator involved, where m = 1 coincides with the 
standard representation of the FL [9] [1(J| [HI 12, 13] . We show by means of a regularisation rule that these 
new representations of the FL are all equivalent to a single unique regularized representation of the FL holding 
for a > 0. We refer to this regularized representation of the operator —(—A) 2 as Fractional Laplacian (FL) 
despite it holds for both, fractional and integer powers of the standard Laplacian. 

The regularization elaborated in the present analysis concerns even (isotropic) functions with respect to 
the hyper-singular point. This regularization makes sense for “functions under the integral”, i.e. we deal with 
distributions (generalized functions) in the generalized sense as presented by G’elfand and Shilov [3]. The 
regularization concerns tempered distributions being defined in the Schwartz space (space of functions which 
have a Fourier transform, see e.g. [16]). 


1.1 Standard representation of Fractional Laplacian 

First of all let us evoke the standard representation of the FL —(—A) 2 of the multidimensional space IR n 
defined for sufficiently good functions <F as (e.g. mmmmmmm) 


- (-A)f$(x) 



<L(x + r) + <f*(x - r) - 2<f*(x) 

r n+a ^ r 


0 < a < 2,n G 111 (1) 


where r • r = r 2 and C( n Q ) > 0 is a positive normalisation factor. With “sufficiently good functions” we mean 
smooth fields <f> with the properties 


lim |<F(x + r) + <F(x — r) — 2<3?(x)| ~ r 2 -> 0 


lim 

r—> 00 


= 0 


(2) 


2 which are in a physical discrete spring model picture equivalent to a self-similar distribution of interparticle springs decaying 
with interparticle distance. Allthough a = 0 is due to translational invariance physically a forbidden value, the regularized 
representation of the FL deduced in this paper holds for all a > 0, (including a = 0). 
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This relation means that a sufficiently good field $ fulfills subsequently introduced general conditions ( |36| ) 
(where a = 0<a<6 = 2) defining the admissible function space in which the standard representation of the 
FL ([Tj) is well defined. We observe that with ([2]) the r-integration of [I]) remains finite at infinity (r = |r|), so 
that 


r const 


(3) 


where we plugged in Qi, and at infinity the r-integration of 0 behaves as 


/ const 


<h(x + r) + <h(x - r) - 2<f>(x) 


dr | rs_/ | — 2^)(x)r 


—a |oo 

const I 


< oo if a > 0 


(4) 


where we plugged in ([2])2- It follows from (J3]) , that the standard FL ([!]) exists in the interval 0 < a < 2 
for sufficiently regular fields <J>, i.e. which fulfill (pi). 

We follow in this paper the convention that the FL is formally defined as operator with Fourier represen¬ 
tation — k a = — (k 2 ) 2 < 0 (k = k|). The sign convention of the FL is here defined in such a way that it has 
all good properties of a Laplaciaifj We come back to the crucial properties of the FL in the course of this 
paper. The normalization factor is obtained as, e.g. 1261 Q3| 


'(n,a) 


2 a ~ 1 aT( c ^) 


7T 2 


r(i-f) 


>o, 


where the F-function is defined by [2] 


0 < a < 2 , 


(5) 


a! 


— r(a + 1) — 


r a e~ T 


dr, 


a > — 1 


( 6 ) 


which exists for a > — 1 and takes for a E IKlo the values of the usual integer factorial with 0! = 1 which 
follows from its recursive behavior al = a (a — 1)!. 

The integration in (l]) is performed over the entire multidimensional space IR n . Another useful equivalent 
representation of the FL is 


- (-A)fcD(x) =: C (n , a) J 5^-^Md-r, 0 < a < 2,n E IN (7) 

where n E IN denotes the dimension of the physical space. In this paper we refer to ([Tj) and ([7]) as the standard 
representation of the FL which exists for sufficiently good functions for 0 < a < 2. 

2 Deduction of the FL from a self-similar linear spring model 

In order to motivate the point of departure physically, we briefly evoke some principal aspects of the self¬ 
similar spring model introduced in [9J. To this end consider a spatially homogeneous distribution of mass with 
constant density po- Each mass-point is represented by the spatial coordinate x and is only connected with 
mass points of distances l s = ha s ( h,a > 0) located at x ± ha s where s £ Zo passes all positive and negative 
integers. We note that by construction the entire set of connecting distances (s E Zo) is invariant under 
the rescalings { l s } —> {a p l s } with p E Zo being integer powers of a prescribed scaling factor a E fR. Note 
that a and a -1 define the same set of sequences { l s } since only s —> —s is exchanged. The harmonic springs 
connecting the material point x with the material points x ± ha s are assumed to have the spring constants 
foa~ Ss = fo{l s /h)~ s , i.e. scaling with a power law with respect to the distances l s . It appears necessary that 
only 5 > 0 is physically admissible , as the spring constants should decay with increasing particle distance l s 
and vanish for ( s —> oo (s -» oo). 

Let u(x) denote the displacement field associated with material point x (or in a generalized context, 
the field variable associated with x). Due to the continuous spatial mass distribution the described harmonic 
spring system can be conceived as a hybrid between lattice and continuum: With respect to the spatial physical 
properties it is a continuum, i.e. all physical quantities depend on the continuous variable x, however, as we 

3 In some references the FL is defined by an opposite sign. 
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will see, the distribution of springs versus distance is a discrete self-similar distribution. The elastic energy 
density associated to mass point x is proposed in the form of a Weierstrass-Mandelbrot type function [9j 


W(x,h) = j ^2 a 5s |(w(x + ha s ) — u(x)) 2 + (u(x — ha s ) — n(x)) 2 | (8) 

s =—oo 

which converges for sufficiently good fields u(x) in the interval 0 < 6 < 2. The total elastic energy is given by 

/ OO 

W(x,h)dx (9) 

-OO 


Let us put here for our convenience the mass density po = 1 and spring constant fo = 1 being equal to 
unity. The prefactor of 1/4 takes into account the double counting of interparticle springs. ([§]) is invariant 
by replacing a —>• a -1 , therefore we can confine without loss of generality to a £ IR > 1. The elastic energy 
density ^ is constructed such that it constitutes a self-similar function with respect to h, and fulfills the 
following self-similarity condition [9] 


W(x, ah) = a°W(x, h) 


( 10 ) 


which implies validity of this relation for all integer powers a s (s € Zo). Self-similarity of a function A (ah) = 
a s A(h) is equivalent with Log-periodicity (ln(a)-periodicity) of the function g(ln(h)) = \n(h~ s A(h)) = g(ln(h) + 
ln(a)) with respect to In (h). Log-periodicity is a typical footprint of discrete scale-invariances observed in 
turbulence. An interesting speculative study on this subject and its link to the FL has been presented recently 
by Chen [20] . 

A function A (h) which is self-similar in the sense (10) is refered to as self-similar at the point h = 0. The 


corresponding notion of self-similarity at a point probably was coined by Peitgen, Jurgens and Saupe [6] . We 
notice that the self-similar elastic energy density ([8]) can be in some cases a fractal function with respect to 
h. The self-similarity condition (10) indicates the absence of a characteristic intrinsic interaction length scale. 
The variable h is a dimensional parameter having dimension of length, however, without the meaning of a 
characteristic (intrinsic) length. 

We define now a ’’Laplacian” by application of Hamilton’s variational principle on ([9]) for the field u which 
leads to 


6V 


Hs , h) u(x)-: ^ 


( 11 ) 


where j- denotes the functional derivative with respect to the field u. We obtain in this way the self-similar 
Laplacian 


oo 

^■{S,a,h) u ( x ) = a ~ Ss ( u{x + ha s ) + u(x — ha s ) — 2 u{x)) , 0 < 5 < 2 (12) 

s=— OO 

converging for the same sufficiently good fields u(x) as ([8]) in the interval 0 < <5 < 2. We refer to this operator 
as self-similar Laplacian since it produces a self-similar function with the scaling property (self-similarity 
condition) 


^(<5, a,ah) ® S,a,h ) (13) 

holding only for a given a and its (positive and negative) integer powers a s s £ Zo- This relation means 
nothing else but that changing of the scale h' = a n h does not change, it just rescales the variables: The 
equation of motion is obtained from Hamilton’s principle by relation © and takes the form 


dt 2 


u(x,t) = A^ a}h) u(x,t) 


(14) 


We refer this equation to as self-similar wave equation. The dispersion relation is determined by A^ a h ^e lkx = 
—1,jJ (5 a)(k) elkx hr the form [9] 


5 





0 < 6 < 2 


(15) 


^.o)^) = 4 H 


a~ Ss sin 2 


kha s 
V o /’ 


This dispersion relation is a Weierstrass-Mandelbrot function (examples are plotted in the figures) con¬ 
verging also in the range 0 < 5 < 2. In the interval 0 < 5 < 1 this function has a fractal graph of Hausdorff 
dimension D = 2 — 5 > 1 which is nowhere differentiable In the range 1 < 5 < 2 ( 15) is a regular non¬ 
fractal (once differentiable) function. The fractal erratic behavior of (15) appears only in the range 0 < 6 < 1 
and is caused by the critically slow decay of the powers a~ Ss as s —> oo. For small exponents 5 —> 0 high 
frequency oscillations are only little suppressed in the series (15) leading to more and more erratic features 
with an increasing fractal Hausdorff dimension D = 2 — 5 which approaches the plane filling dimension 2. In 
its entire range of convergence 0 < 6 < 2 the Weierstrass-Mandelbrot function type dispersion relation (15) is 


self-similar with respect to the non-dimensional wave number kh , namely a ^(akh) = a°wig a \{kh). 

In the two figures two dispersion curves are plotted, for fractal regime (erratic curve for first figure) 
0 < 5 < 1 and non-fractal regime 1 < 5 < 2, respectively (erratic curve of second figure). The smaller 5 
the more erratic become the dispersion curves. The smooth curves indicate their fractional continuum limits 
a —> 1 (given in explicit form in below equation (23)) where the erratic fractal characteristics are wiped out. 


2.1 Fractional continuum limit 

An important limiting case which we refer to as the fractional continuum, limit is obtained for a —5 1 in 
relations ([8]), (12) and (15) [13!]. This limiting case connects the self-similar discrete spring model of the 
subsequent section with fractional calculus, especially the discrete ID self-silmilar Laplacian (12) with its ID 
FL counterpart. To demonstrate this in brief, let us define the fractional continuum limit by means of a 
self-similar function A a (h) of Weierstrass-Mandelbrot type 


A a(h)= £ 


-,—Ss 


f(a s h) 


(16) 


which is self-similar with respect to h fulfilling A (ah) = a s A(h) and converges for sufficiently good functions 
/. We define the fractional continuum limit by 


lim | ln(a)|A a (/i) = h° 

a —>1 


/( r ) 


/ 0 T' 


<5+1 


dr 


(17) 


which remains finite for o—>-l. We notice then the following observation: The physics described by the 
FL defined by the fractional continuum limit (17) of the self-similar Laplacian is up to a prefactor ~ h s 
independent on h. The constant h has therefore only the meaning of a dimensional normalization factor 
(having physical dimension of length), but does not have the meaning of a characteristic length-scale which 
is nonexistent for self-similar functions fulfilling A (ah) = a s A(h). We can perform the fractional continuum 
limit for any self-similar function of the form ( |16[ ) and the range of convergence for 5 is the same for both, the 
discrete case (16) and the fractional continuum limiting case ©• 

Further illuminating is to consider the fractional continuum limit of a self-similar density function p a ^ 
generated by Dirac’s d-functions 


Pa,h( T ) = \ a SS ^ T ~ haS } f = 

0-51+ ' 


h S 

£t (S+1 


(18) 


where in this limiting case 0 < £ = |ln(a)| << 1 is assumed to be a positive “small” but finite quantity 
(a ~ 1 + C). The fractional continuum limit of (16) can then be represented by 


roo 

Aa(/l) = / Pa,h( T )f( T ) d7 
JO 


(19) 


assuming in the fractional continuum limit the representation ©• Performing the fractional continuum limit 
for the self-similar Laplacian (12) defined by 
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( 20 ) 


lim A { 5 , a ,h) = A m 


we obtain (up to a positive prefactor) the standard representation of the FL in one dimension [10, 11, [12] 

h s u(x + r) + u(x — t) — 2 u(x) 


= 2Q L 


15+1 


-dr, 


0 < 6 < 2 


( 21 ) 


where the additional prefactor \ comes into play by choosing the integration limits ±oo and by taking into 
account that the integrand (21) is an even function with respect to r. The standard representation of the 
FL (21) exists in exactly the same interval 0 < 6 < 2 as its discrete self-similar counterpart (12). It is now 


relatively straight-forward to generalize (21) to the multidimensional case of n dimensions of the physical 
space 


A 


(a,C)«( x ) = ^ 


h a f -u(x + r) + rt(x — r) — 2it(x) 




d n i 


0 < a < 2, n G 


( 22 ) 


where x£ R" and the integation is performed over the entire IR”. The generalization (22) to the IR n includes 
(up to the prefactors) the standard representation of the FL ([Tj) . The generalization (22) of (21) maintains the 
spatial isotropy and the dimensional prefactor h a in (22) maintains the physical dimension of the self-similar 
Laplaciarj^] (dimension of the energetically conjugate quantity of the held variable) independent of dimension 
n. The dispersion relation (15) takes in the fractional continuum limit the form of a power law m 


Jj{k) = A 5 k 6 


(23) 


which is in the same time the dispersion relation obtained from the continuum representation (21) where the 
constant As is obtained as m 


As = 


h 4 5 


7T 


C r(<5 + 1) sin f 


> 0 


(24) 


where the relation to the normalization factor (5) (dimension n = 1) is A a = —. It is important to note 

that representation (22) exists only for exponents 0 < a < 2 for sufficiently regular fields tt(x). We emphasize 
that this restriction is a property of the spatial standard representation only and not of the operator of FL 
itself which should exist with all good properties of a true Laplacian for any positive a > 0. 


2.2 The generalized self-similar potential - towards a generalization of the FL 

In this section we introduce a generalization of the above self-similar spring model. Hamilton’s variational 
principle defines then in a rigorous a manner a self-similar Laplacian which takes in the fractional continuum 
limit for fractional non-integer powers new representations of the FL. To this end it is convenient to introduce 
the shift operator by 


D x (h)u(x) = u(x + h) (25) 

where we consider first of all the ID case. For sufficiently smooth (infinitely differentiable) fields u(x) where 
it is sufficient to consider only scalar fields u, the shift operator has the representation 

a 00 h m d m 

O l( ft ) = e x p(A-)= — (26) 

m=() 

which generates the Taylor serie^] We notice in brief some essential properties: D(h = 0) = 1 is the unity 
operator of zero-shift and D^(h) = D(ph) for any ju G IR and the shift operator is unitary, i.e. D(—h ) = 

4 whereas the physical dimension of the FL is always 1/ (length) 01 

5 We notice that since the shift operator is an operator function of its generator 4-, it has the same eigenfunctions e lkx where 
D x (h)e ikx =e ikh e ikx . 
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D l (h) = D\.(K) where with (,.)t we denote the adjoint operator of The shift operation fulfills all 
properties of a continuous commutative Abelian group. 

Let us now introduce a generalized form of discrete self-similar elastic potential ([8]). The generalization is 
such that it involves not only two-point but self-similar harmonic m + 1-particle interactions which come into 
play by quadratic forms of finite differences of orders m = 1,2,.. e 111, namely 


r OO 

Wm(x, h) = E a ~ dS {[ D x (ha s ) - 1 ] m u(x)} 2 , 


0 < 5 < 2m 


(27) 


Note that this series is a function of type Weierstrass-Mandelbrot where self-similarity is reflected by 
W m (x,ah ) = a s W m (x,h). This self-similar series converges for sufficiently regular fields u(x) in the interval 
0 < 5 < 2m. For m = 1 (27) coincides with the potential ([8]), i.e. with the self-similar chain model previously 
introduced [9]. Let us again put f m = 1 subsequently. The total elastic energy functional is then given by 

/ OO 

W(x,h)dx (28) 

-OO 


self-similar Laplacian is obtained by Hamilton’s principle 


The generalized self-similar potential (27) fulfills also the self-similarity condition (10) and a generalized 

(29) 


A n- sv ™ 

A {mAh) u(x) . §u ^ 

Then this self-similar Laplacian is given by the representation 

OO 

A( m ,s,h) u ( x ) = E a~ 5s A 2m (a s h)u(x), 

s=— OO 

where the finite difference operator of even order 2m comes into play 


0 < 5 < 2m 


A 2m(h) = (-1 ) m+ \D x {±) - D x (~)) 2m = - (2 - D x (h) - D x (-h)Y 


m e m 


(30) 


(31) 


which is symmetric A 2 m (h) = A 2 m (—/i) and hence self-adjoint. We restrict ourselves in (27) on finite differ¬ 


ences of integer order m 6 111. However, we emphasize that the entire analysis to follow remains true and can 
be performed in full explicit form also for finite differences of fractional orders, where the orders m must be 
strictly positive m > 0. The positiveness is the only real condition which is necessary to maintain translational 
invariance (eigenvalue zero at k = 0) for the FL to be deduced. The derivation of ( |30[ ) involving operator © 
is given in details a recent paper in the context of non-local constitutive models |14j. Some important but 
elementary features of the difference operator which are used to obtain (|31|) are outlined there. 


We further notice that for the case m = 1 the self-similar Laplacian Q12[) is re-obtained. (30) converges 
in the same interval 0 < S < 2m as the elastic energy density (27). Taking into account that e lkx are 


eigenfunctions of the shift operator and hence of (30) we obtain the dispersion relation from A ( mi s,h) elkx = 
( kh)e lkx in the form of Weierstrass-Mandelbrot type functions 


—w; 


m,S V 


kha s 


<i(^) = 4 m E a ~ SS sin 2m ( - ), 


0 < 5 < 2m 


(32) 


involving powers of orders 2m of the sine. The self-similar Laplacian (30) as well as the elastic potential 


(27) and the dispersion relation (32) are self-similar functions in the sense of a condition as (10). For m = 1 


this model coincides with the self-similar spring model of the previous section leading to Laplacian (12) 


The dispersion relation (32) converges as well in the entire range 0 < <5 < 2m and is a function of the 


type Weierstrass Mandelbrot which is nowhere differentiable and with fractal features within 0 < 5 < 1 (see 
the figures representing the case m = 1). The further analysis of the fractal characteristics of Weierstrass- 
Mandelbrot type functions, however, is not subject of this paper. For further aspects of those functions we 
refer e.g. to jjT]. 

















The fractional continuum limit a 
then yields 


1 + C (0 < C << 1) as introduced in section 2.1 of Laplacian (30) 


A, 


= h L 


A 2 m(r)u(x) 


15+1 


dr, 


0 < 6 < 2 m 


(33) 


which has the form 


coinciding with (|21|) for rri = 1. For the analysis to follow we will also need its n-dimensional counterpart 

h s f A 2m (r )u(x) 

A m ,s,c u ( x ) = ^J —' 


-d n r, 


0 < <5 < 2m 


(34) 


which is well defined in the entire interval 0 < 5 < 2m where the symmetric difference operator A 2rn (r) is 
defined below (eq. ©)• Note that ( |33[ ) and ( |34| ) are non-singulai integrals which exist within the entire 


interval 0 < <5 < 2m. The hypersingularity at t = 0 is removed by the tempering property of the difference 


operator (31) behaving as ~ r 2m for r —> 0. We will show in the subsequent analysis that (33) and (34) 


represent up to normalization constant the respective ID and nD FLs for fractional exponents 0 < a < 2m. 

However, for integer exponents ^ E HI the non-regularized expressions (33) and (34) represent integer 
powers of the Laplacian A only for a restricted subspace of tempered functions and looses otherwise its 
analyticity by the appearance of not well defined integrals. This analytic anomaly at integers merits a deeper 
analysis which is beyond the scope of this paper. 

Subsequently we deduce a regularized representation of the FL which avoids this anomaly representing 
for the full range f > 0 the powers of the standard Laplacian A including the integer powers | G IHq. 


3 New representations for the FL 


Let us consider an integral of the form 


InAQ — 


/(k-r 

y>a.-\-n 


-d n r = I n ,a{k = 1 )k c 


(35) 


which is performed over the entire IR' 1 and we notice that it scales as ~ k a for sufficiently good functions / 
for which this integral exists where the non-zero I n ,a(k = 1) appears as a normalization constant. Let us 
hence assume that /(£) = /(—£) is an even function. The integral (35) is regular for functions / of which the 
hyper singularity 


a— 1 


at r = 0 is removed. Such functions require the following asymptotic properties: 


1/(01 < M b , 

1/(01 < foot*, 


£ 


(36) 


oo 


where the integral (35) is regular (non-singular at r = 0) only if a < a < b and a < b. The first condition 
(36)i has in the integral tempering effect at r = 0 and renders (35) locally integrable at r = 0. The second 


condition (36)2 guarantees integrability of (35) at infinity. The conditions (36) define the admissible function 


space in which (35) is defined. In (36) /o, /oo > 0 denote characteristic constants of the function / and where 


a < b is a necessary condition for / to be admissible. Now, generate a function 

/ 2 (h) = A 2 (h)ti(x) 


(37) 


where A 2 (h) indicates the symmetric second order finite difference operator of (31) for m = 1, then we notice 
that / 2 (h) fulfills with respect to its dependence on h the properties (36) with a = 0 and 6 = 2. Since the 


standard representation ([Tj) of the FL is constructed by a second difference operator of the form (37), it indeed 
is regular within the interval 0 < a < 2. 

As a generalization we can employ a function / 2m = A 2m (r)u(x) generated by the symmetric finite 
difference operator of order 2m and define the integral 


A 2m (r)u(x) 


r n+a 


d n r 


(38) 
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which is regular for a = 0<a = b = 2m. Further it follows from (35) that this integral has for it(x) = e* kx 
the scaling behavior of (35) 


— T mna (k — 1)^ 


0 < a < 2m 


(39) 


The Fractional Laplacian (FL) —(—A) 2 is defined as the operator which has in Fourier space the (negative 
definite) representation — k a . We will demonstrate that we can use integrals of the form (38) to derive spatial 
IR^-representations for the FL which represent powers of the standard Laplacian A. We emphasize that the FL 
is a uniquely defined operator, so that all spatial representations which generate the FL need to be equivalent 
and coincide. 

We propose hence the following representation for the FL which will be later justified by its exact definition 
by Hamilton’s variational principle 


m ( A) 2 li(x) — C)ra,n,a J 


*2771 


»u(x) 


r a+n 


d n r, 


0 < a < 2m, 


(40) 


which exists for 0 < a < 2m for sufficiently regular function^ it(x). In (40) we have employed the symmetric 
(self-adjoint) finite difference operator of order 2m (m > 0), namely 


i2m 


f r r 'l 2 tti 

(r) = _ {2 -D x (r)-D x (-r)r = (-ir+ I (D(-)-D(--)j , m e N(m > 0) 


(41) 


where D x indicates the unitary shift operator which shifts the vector argument x —> x + r, i.e. D x (r)u(x) = 
u(x + r). The symmetric finite difference operator (41) renders (40) locally integrable at r = 0 due to its 
tempering behavior ~ r 2m for r “small”. The tempering effect is sufficient if a < 2m and is present for 
sufficiently smooth fields. Further we deal with “good fields” which are vanishing at infinity and as a conse¬ 
quence | A 2 m(r)u(x)| ~ -m!m! u ( x ) = cons t converges versus a non-zero constant (i.e. independent on r) for 
a > 0. From these considerations it follows that (40) is a non-singular integral existing for 0 < a < 2m. The 


normalization constant C mna is introduced that (40) has Fourier transform —k a . We show in the appendix 
ra that this normalization factor is well defined in the entire range 0 < a < 2m. However, we will see that 
it is sufficient to evaluate it explicitly for fractional powers ^ ^ IUq only. 


We prove subsequently the following features: Expression (40) represents powers of the conventional Lapla¬ 
cian A only for fractional (non-integer) exponents 9) ^ IKI. However, in the regularized representations of the 
FL to be deduced, integer powers f € I'll are also included, including the (due to translational invariance) 
physically inadmissible value a = 0. The regularized FL gives a spatial representation in the distributional 
sense for the powers of the standard Laplacian — (—A)a (Equation (55) below) which holds for any a > 0. 


Especially this regularized representation for the FL is independent of m and coincides with (40) in the in¬ 


terval 0 < ^ ^ IKI < m. This regularized representation for the FL takes for integer exponents j G Wo a 
distributional representation for the integer powers of the conventional Laplacian A including for a = 0 the 
negative unity operator. 


We will show now that the FL (40) can be reduced to a unique simple universal representation which is 


independent of order m. This reduction requires the introduction of a regularization to allow certain manip¬ 
ulations. This regularization leads to a simple representation of the FL holding for any a > 0 representing 
powers of the standard Laplacian A. To this end we consider the more explicit evaluation of the FL (40) 


777 ( ^) 2 «(X) - Ct77,?7 ,a j 


1 


(2m) 


r a+n j uiluil 


u X 


) + £ 


(2m)! 


—( (m + p)\(m — p)\ 


which exists for exponents 0 < a < 2m. The prefactor C m 
that the eigenvalues of the FL are given by —k a and is deduced in the appendix A.l and is given by 


(—l) p [u(x + pr) + u(x — pr)] ^ d n r 

(42) 

is a normalization constant which guarantees 


6 i.e. infinitely differentiable and vanishing at infinity. 
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(43) 


A r , 


= C. 


-l 


— Un.aVm.a ^ 0 


(- 1 ) 


m+1 7r^r(^) 


2 ar(«±H)r(a+ l)sin^f 
where the integral over the unit-sphere 


(D x (l)-D x (-1)) 


2m | 


a 


|a=o, 0 < a < 2m, - £ 


Ut. . rv — 


|ni| a dfl(n) = 


71—1 , l 1 \ 

27r^r(^±i) 


> 0 a > — 1 


(44) 


'H=i- ' r ( 2 ) 

comes into play and where n\ denotes any Cartesian component of the unit-vector n and with the integral 

2 „,-„ r™ 2ro © J5>0 

(45) 


V — 2 

v 771,.CM - 46 


= (- 1 ) 


m+1 


£oH-l 
7 r 


2 Q+1 a! sin 


d£ > 0 

(D X (1) - D x (-l)) 2m |ArU =0 , 


a , 

0< 2* 


< m 


Note that relation (45)2 holds only for non-integers 0 < “ ^ IN < m. 

It is crucially important to note that (42) does not allow commutation of the binomial sum and r- 
integration, since the r-integrals become singular without the tempering effect at r = 0 of the A 2 m - difference 
operate^ In order to allow this commutation integral (42) has to be regularized. We observe that the in¬ 
tegrand is an even function in r. By introducing polar coordinates r = n?’ (n denotes the unit-vector in 
r-direction) (42) can be further written as 


—m(—A) 2 u(x) = —C m n a [ dft(n) ( I 

J |n|=i \ rn\m\ Jo 


“ dr!®+2 f 

p 


(2m)! 


r*Oi-\~ 1 


(m + p) ! (m — p )! 


(-1) P / 

Jo 


00 , u(x + prn) 
dr-— 


_ (46) 

We emphazize that by commuting integration and summation in (46), the termpering effect of the A 2 m 


difference operator does not take place, and hence (46) is a sum of singular integrals at r = 0 requiring 


regularization in order to be defined. The regularization needs to eliminate the divergence of any integral 
term of the form 

l/^ M ( x + P r )l = N“l/^wf u ( x + r )l^ 00 > a>0 ( 47 ) 


which diverges due to the hypersingularity of order ~ r 
a regularization with the following properties: 


—a— 1 


of the integrand at r = 0. To define (46) we need 


(i) The regularization of U^ ) after the commutation coincides for good fields u with the unregularized well 
defined representation 1 4 A 5 


(ii) The regularization renders (42) finite in order to allow commutation of integration and summation in (46) 


We employ the following regularization rule that fulfills criteria (i), and (ii) which was already employed 
earlier EQ CEDi HU Q21 dg 


-lei 


—a— 1 


sin — = lim K(e — i£) 
2 e->o+ v 


— a— 1 


(48) 


where 3?(..) denotes the real-part of (..). This relation is an exact equality for non-zero |£| > 0. For £ = 0 
the right hand side guarantees local integrability at £ = 0 by removing the hypersingularity of the left hand 
side which occurs when a > 0. We can say that regularization rule (48) is an equation in the distributional 


7 A 2m (r) 


’ r 2m as r —> 0) 
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sense, i.e. it holds “under the integral”. The regularizing effect of (48) can be demonstrated for instance by 
considering the integral (£o > 0) and a > 0 


IregiCo, «) = £ 1™_ ^ 


K o 


( e _ = 


sin 


a 


C-a 

*0 


(49) 


Integrating the left hand side of (48) yields a singular integral with a diverging term-“ (£ ~t 0) 


at the lower integration limit £ = 0. The crucial feature of regularization rule (48) is that the divergence of 
the lower limit is removed leading to the finite value (49). Furthermore (48) gives for integer f G IFI o to the 


left hand side a distributional meaning by means even derivatives of Dirac s 5-function. That is, by employing 
in rule (48) integrals containing — |C| _a_1 sin ™ are well defined. We will show subsequently that 
the FLJ46I) contains such term and hence the FL (46) is well defined also for integers j € N o if regularization 
rule 


is applied and by accounting its distributional form for these cases 


sin 


= lim 


„-a+l 


|£| a +i e^o+ 1 (£ + ie) a+1 


} = (-l) p 


7r d 2 P 

(2 p)\ df 2v 


a 


7T =P G PUf 


(50) 


where 6(£) denotes Dirac’s 5-function. It follows that in the limiting case of integer exponents (50) takes 


completely local distributional representations of even-oder derivatives. Further we notice that in the limiting 
case a —> 0+ the regularized integral (49) tends to a finite value I reg —> %. 


We hence can write for (46) which is defined only by employing regularization rule (48) 


- m (-A)2ii( x ) = 


jm-i dn(n > i w r ^ + 2 1 ; 


u(x + £n 




(2m)! 


£“ +1 “ (m + p)\ (m — p )! 


(-1 Yp° 


(51) 


In this way we can for the terms p / 0 separately re-scale the integration variables to £ = pr and obtain 
for all terms of p ^ 0 the same integral with integrand ■ The term for p = 0 contributes only as 


integral of . This term p = 0 is vanishing in the regularization (leading to an integral of the form (49) 
with £o —>• oo). So the constant term ~ u(x) which corresponds to p = 0 does not contribute. Hence we can 
rewrite the sum which can be extracted from relation (51) a: 


^E 


(2m)! 


( (m + p)\(m — p)\ 
(- 1 ) 


(-i)V = §Wo“ + 2]T 


(2m)! 


m!m! 


—( (m + p)\{m — p)\ 


(-1 ) p P° 


(52) 


[D x (l)-D(-l)] 2 m \\n x=l 


a > 0 


where the vanishing term of p = 0 has been added for convenience. Note that this relation is in the range 
0 < < m non-vanishing only for ^ ^ fkl. 

We obtain then for the FL the following (without regularisation rule (|48|) singular) representation 


— m (-A) 2 «(x) = 

(-ir +1 C m ,„,„2-« [D A {1) - D(-l)] 2m |A|“| a _ 0 J + r) 

f d n r , 

■J-a TT M ( x + r )’ a>0 

u( r) 


2a! sin ^ f d” 


( 53 ) 


7 tU, 


n,a 


-(-A) 2«(x) = C n , c 


_ r \n-\-a 


X — r 


d n r, 


a > 0 


8 with the shift operator Dx(a)f(X) | A =o = /(A + <i)|a=o = /(a). 
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which is independent of m and well defined in the distributional sense when employing regularization rule 
(48). The normalization constant which is indepedent of m is then obtained by 


Cri..rv — 


2a! sin f T(^)a! sin 


irU n 


7T 




(54) 


coinciding for 0 < a < 2 with the normalization factor of the standard representation ([5]). Then we obtain 
the regularized representation for the FL 


-(—A)f u(x) = 
-(—A)2 u(x) = 


— lim Vt(i a+1 f n(X + r) d " 
vrC/ n ,«^o+ V J (r + ie) a + l r"” 1 


hm ^- 1 


7r"2 1 r(^±i) e_s>0+ 


u(r) 


d n r 


(55) 


(|x — r| + ie) a+1 |x — r| 


n— 1 


a > 0 


We emphasize that (55) holds for all positive exponents a > 0 including integers (with zero) that (55) 
represents the powers of the convential Laplacian A. We notice that for fractional powers % INo, the 
regularized FL (55) is a non-local convolutional (scale-free and self-similar) operator which remains non-local 
whatever the scale of r representation may be, whereas (55) becomes a local operator for integer powers 

Mo 

For 


- e 
2 


can be considered as is the analytical extension of the FL (40) to full range of exponents a > 0. 
D case n = 1 the regularized FL of (50), (55) yields 


d 2 
dx 2 


)2 u{x) = lim — — / 

e—>0+ TT [ J- 

° u(x + t) 


u{t) 


a! sin if 

IT 


la+1 


-dr, 


-oo (x - T + ie)" +1 


a 


d t 


a > 0 


(56) 


a > 0,- i M 

d 2 


where integers | = p £ N o yield the required result — (— j^) p - Relation (56)2 is its unregularized singular 
representation which is defined only in the sense of (56)i, i.e. by employing regularization (48). 


Let us now analyze in full details the case of n dimensions for j = p G IN o to prove that for n dimensions 
the regularized representation of the FL (55) indeed yields for j = p G INo the integer powers — (—A) p of the 
conventional Laplacian. To this end we can write for (55) by using the shift operator (26) where only the 
even powers contribut^] 


- (-A)2 =p u (x) = - 


1 


E 


(n • V x ) 2s u 


Un,2p o—n |n|=! 


(2s)! 


( TT Jo 


„2 s 


o (r + ie) 2 P +1 


dr 


(57) 


where U n ^ P is the surface integral (44) for even integers a = 2p £ INo- Further we observe by taking into 
account relation (56) that the r-integral yields after 2 p partial integrations 


2 


(2 p)\i 


|,’2p+l r oo 


„2s 


tt Jo (r- + ie) 2 P +1 


dr 


(58 > 

roc 

= 2(-l) p (2p)\5 ps / 6(r)dr = (-l) p (2 p)\6 ps 
Jo 

being non-zero only for the order 2s = 2 p ( 5 ps denotes the Kronecker-<5). With this relation we can write for 

(5^ 

u The gradient V x acts only on the x-dependence of u(x). 
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(59) 


-(-A )*u(x) = 


(_1)P+! 


u. 


n,2p |n|=l 


(n • V x ) 2p dfl(n) u[ 


x 


To evaluate this surface integral, we utilize (44) to arrive at 


J| n |=i dft(n)(a • n) 2p = U n , a=2 p (a • a) p 


— Un,a=2p Y. 


(2 p)\ 


J pi n 2 s\.yp 


(60) 


Vi+P2+-Pn=P (2pi)!(2p 2 )!..(2p„)! 1 2 


where a denotes a constant vector (independent on the surface unit-vector n). In (60) only even combinations 
of powers of the Cartesian components n* of the unit vector yield non-vanishing contributions. Replacing in 
this relation the monomials by partial derivatives a\ Pl a^ 12 ..a%P n = 2pi ^ - 2pn gives for (59) finally 

- (-Afu(x) = (-1)P +1 (V- V) 2p , pgM 0 (61) 


which indeed is the correct integer powers of the standard Laplacian including for the case p = 0 where © 
yields — u(x) since (—A) p 0 = 1 yields the unity operator. The sign (—l) p+1 alternating with order p appears 
always such that (61) has Fourier representation — k 2p with all good properties of a Laplacian. 


3.1 Hamilton’s variational principle: the FL versus self-similar Laplacian 

In this section we demonstrate that harmonic self-similar elastic energy density functionals of Weierstrass- 
Mandelbrot function type lead via application of the fractional continuum limit defined in section |2.1| tog ether 
with Hamilton’s variational principle to rescaled versions of the non-regularized FL of the type (|40|) with 
limited intervals of existence 0 < a < 2m. Application of regularization rule (48) then allows to arrive at the 


universal regularized fractional Laplacian FL (55). We show this to be true for any admissible self-similar 
harmonic potential energy densities of Weierstrass-Mandelbrot type. In this way we demonstrate the general 
interlink between self-similar elastic potentials (describing physically the absence of an intrinsic lengthscale of 
the range of interparticle interactions), and the appearance of the regularized FL (55) as an universal footprint 
of self-similar elastic systems. Our goal is to demonstrate in this section that the physical information on 
the details of the harmonic interactions of the self-similar potentials appears only as a characteristic scaling 
factor in the rescaled FL. In other words: The physics on such self-similar harmonic systems is independent 
on the details of the self-similar potential. This characteristic scaling factor does not contain any absolute 
length scale being again a consequence of self-similarity (absence of internal length scale of the interparticle 
interactions). 


Let us consider a discrete elastic energy density of self-similar Weierstrass-Mandelbrot function type (27) 
which can be written as 


1 


W a (x,h) = - ^2 a as {(D x (ha s ) - l) m u(x)} 2 , 


0 < a < 2m 


(62) 


which fulfills the self-similarity condition W(ah) = a a yV(h). Application of Hamilton’s variational principle to 
(62) yields the self-similar discrete Laplacian (|30|) which is also of Weierstrass-Mandelbrot type. The fractional 


continuum limit a —> 1 of (62) yields by using Ql7|) the ID integral (a=l + C>0<C~ ln{a) « 1) 

h a r°° dr 


roo 

W a=1K (x, h) = -J o ^ {( D x (t ) - 1 ) m u{x)Y 


2( Jo r a+1 

The total elastic energy is then given by 


0 < a < 2m 


(63) 


V r 


m,l,cx. — 


Wa=i+t(x,h)dx 


— . - 


h a r°° r°° drdx 


(64) 


Af I I q~' 

J —oc J — oo ' 


q ;+1 


{(D x (r) - l) m u(x)} 2 
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where the additional prefactor i comes into play by taking both integrals in the limits (— 00 , 00 ). The 


generalization of the ID self-similar elastic energy (64) to the n-dimensional space yields the representation 


h a 


{(Ac(r) - 

y*a-\-n 


d n rd n x, 


0 < a < 2m 


(65) 


Application of Hamilton’s variational principle to (65) yields the self-similar continuous Laplacian (|34l) which is 


up to its different prefactor coinciding with the non-regularized FL (40). The positiveness of this normalization 


factor (43) appearing in (40) reflects elastic stability (positive definiteness of the self-similar potential). 


Performing the procedure of the last section toegther with regularization rule (48) finally yields a rescaled 
version of the regularized FL (55). 


It is noteworthy that, the normalization constant (54) of the regularized FL is due to the oscillations of 


sin Tr alternating in sign which indicates that the regularized representation of the FL (|55|) cannot be derived 


directly from a variational principle as a functional derivative of a potential such as (61 

Now let us analyze a generalization of the elastic energy functional (65) in the following fornf^j 


Va[u] = i 


^77. j, _^ 

—^r E v \P-q\< y + p T >(y + <? r ) ^ °’ 

p,q 


0 < a < f3 


( 66 ) 


where we assume that (J 66 | is well defined for 0 < a < (3 where (3 depends on the choice of the matrix V\ p _ q \. 
Further we have V a [u = const] = 0 which indicates translational invariance (zero elastic energy for uniform 
translations u = consrH and V a [u] >0 otherwise). The self-similar continuous “Laplacian” Ay of this system 
is defined by Hamilton’s variational principle leading to 12 ' 


Ayti(x) = — 


<7u(x 


■V a [u] 


1 f d n r _^ 

Ayu(x) = - - / —^ Y, V pq (u(x + (p - q)r)) + u(x - (p - g)r)) 
J p,q 


(67) 


AWx) = -- 


d n r 

p n+a 


E V pq U (* + (P “ ^) r ) 


P>9 


To evaluate last step (67)3 we face the same divergence problem as in (47) when we commute summation 
and integration. To render this commutation admissible we need to apply regularization rule (48) which 
regularizes each singular integral term of the sum. Then (67)3 takes the representation 


1 


Ayu(x) = ~^Y V l 


pq 


d n r 

r n+a 


u(x + (p — q)r) 


( 68 ) 


pq 


which is as relation (56)2 a sum of singular integrals, and only defined by employing regularization rule (48) 
By plugging in (48) into ( 68 ) we can integrate any term of the sum by rescaling its radial integration variab 


(by introducing for each term p 7 ^ q the new intgration variable \p — q\r = £ —> r) and with the same reasoning 
which lead us to to above relation we can evaluate each term of the sum ([ 68 ]) by using 


tt(x + (p— q) r) 


d n r = \p — q\ a J 


u(x + r) 


r n-\-a 


d n r 


(69) 


J /pTl-\-(y. J rp l I 

This relation is for a > 0 also valid for p = q, since the regularized value of the integrals occuring at p = q 
are vanishing due to (49) 

f d n r 

/^TT=0 (70) 


1() Where the symmetric constant matrix V pq = V qp = Vj p _ 9 | having only one single eigenvalue zero, due translational invariance 
with V(u = const ) = 0. 


11 where = const is not a “good field” in the sense of (361 

^"where we write here just the integrals in their singular representations which are to be defined by applying regularization rule 
(481. 
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So we obtain for the Laplacian (67) just a rescaled version of the FL (53), namely 

u(x + r) 


f 

Ayti(x) = —Ay / 


r n-\-a 


-d n r 


(71) 


being defined only with regularization rule (48) and in this way being a rescaled version of the regularized FL 
( |55| ). We emphasize that ( fTIj ) is only definedm regularized sense, i.e. by plugging in regularization rule (48). 
The characteristic positive scaling factor 


-Av = n V \p-q\\P~ <l\ a > °> 


a > 0 


(72) 




can be conceived as a self-similar stiffness matrix containing the “scale free” information of the matrix V \ p - q \, 
but it does not contain any information on an absolute length-scale of the interparticle interactions. This 
is expressed by the self-similar power functions |p — q\ a : Any change of scale of Jn — q) —> X(p — q) only 
rescales the characteristic scaling factor (72) by |A|“. For the special potential (65) the scaling factor (72) 
is evaluated in the appendix A.l From this demonstration follows that the regularized representation (55) 


of the FL is the universal and unique charactersitic operator of all harmonic self-similar systems with elastic 
potentials of the general form (66). Note that the scaling factor ( |72| ) is only well defined for the physically 
admissible exponents a > 0 which reflects its “physical origine” (respecting translational invariance) whereas 
the regularized FL (55) is mathematically well defined in the distributional sense for a > 0 including also the 


physically forbidden value a = 0. 


4 Conclusions 

The main goal of the present ana 
which exists for any a > 0 (Eq. I 
tion (relation (|48|) ). We analyzec 


ysis was to deduce a spatial distributional representation for the FL — (—A) 2 


55)). The deduction of this reprepresentation inevitably requires regulariza- 


self-similar elastic potentials of the form (66) which include potentials 


that are obtained by the fractional continuum limit from discrete self-similar Weierstrass-Mandelbrot type 
potentials of the for m (|62[ ). Hamilton’s variational principle then defines in a rigorous manner self-similar 
Laplacians © and (|34[), respectively which take the representations of rescaled variants of the universal 
regularized FL (55). Their characteristic scaling factors contain the information on the self-similar elastic 


potentials (stiffness matrices). Since the characteristic information on the elastic potentials appears only in 


the form of different scaling factors of the regularized FL (55), this operator has universal character for har¬ 
monic systems with self-similar, scale-free interactions. All these self-similar harmonic systems are uniquely 


governed by the regularized FL (55). The appearence of this universal regularized representation of the FL 
opens a considerable reduction of difficulty in order to handle physical problems which are described in a frac¬ 
tional calculus context. The extended range of validity to exponents a > 0 especially opens newly emerging 
problems to the analysis which are described by exponents outside the Levy-interval (0, 2) such as recently 
described for distributions of fluctuations in financial market indices [5]. 

The self-similar spring model invoked as point of departure for our approach as well as the generalized dis¬ 
crete self-similar potentials of Weierstrass-Mandelbrot type (62) provide the interlink to the inherent fractal 


properties of these systems, having as a consequence dispersion relations (32) of the type of Weierstrass- 


Mandelbrot functions with fractal features in the regime of exponents 0 < 5 < 1. In this way we have demon¬ 
strated the relation between fractality of certain physical properties such as of the Weierstrass-Mandelbrot 
type dispersion relation (as visualized in the figures) and the “generalized” fractional Laplacian approach 
presented in this paper. There is a big deal of interdisciplinary scaling invariant problems in various physical 
contexts such as in critical phenomena, anomalous diffusion and related stochastic motions such as in tur¬ 
bulence, where the regularized representation of the FL (55) and its links to Weierstrass-Mandelbrot fractal 


functions, may be useful. We hope the present paper inspires further works on those and related problems. 
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A Appendix 


A.l Explicit determination of the normalization constant n Q 


In this appendix we determine the normalization constant which comes into play to normalize the 

representation of the FL (40) such that its Fourier transform is — k a . in explicit form. The normalization 
constant is determined by the exigence that the FL has Fourier the transform —k a . Accounting for this, it 
follows from (39) that 


C 


m,n,a '-7 m,n,a(k — 1) — 1 


For our convenience we put —I m , n ,a{k = 1) = A m . n ^ a = C m 1 na which takes the form 

( e #_ e -#)2™ o2ro _ a f sin 2m (fc • r) 


A — (—'\Y 

^777,77,0! — V -*-/ 


y.n-\-cx 


-d n r = 2 


r n-\-a 


-d n r > 0 , 


0 < a < 2m 


(73) 


(74) 


where k ■ k = 1. The constant A m ^ n ^ a > 0 is positive. Its positiveness is crucial for the FL in order to have 
all good properties of a Laplacian. We notice that the constant (74) is defined for any 0 < a < 2 m, i.e. also 


for | E III. By introducing the polar coordinates r = n (r = |r| and n • n = 1) and choosing the coordinate 


system such that k ■ r = rn\ where ni is a director cosine. We notice that (74) can be factorized according to 
{Am,n,a = C mna ) 


= C. 


-l 


— Ur..ryVrr 


^m,n,a — ^n,a r m,ai 

with the surface integral over the unit sphere 


0 < a < 2m 


(75) 


/■ 2'K'°A^T( g AA\ 

Un,a = |ni| Q dP(n) = — + , 2 > o (76) 

J|n|=l 1 1“2-J 

This integral U n>a exists for a > — 1 and is explicitly evaluated in the subsequent appendix A.2 For 
instance (76) takes in ID the value U n= \ a = 2. In (75) we further have the integral 


V m ,a = 2 


2m—a 


siiL m (0 

£c*+i 


d£> 0 , 


0 < a < 2m 


(77) 


The existence of (77), especially 0 < V m ,o. < 00 and hence A m n ^ a > 0 which indicates “elastic stability 


(positive definiteness) of the system which is due to the fact that (34) is the outcome of Hamilton’s variational 
pronciple of a positive definite elastic potential in the entire open interval 0 < a < 2m. A further point is 
crucial: We see in view of (76) and (77) that the normalization factor as well as (40) is a well defined operator 


for any 0 < a < 2m, especially also at integer points |=pElll. The normalization factor C m ^ n ^ a of (40) is 
given by (75) 


Am,n,a=2p — C m n ,2p — U n ,2pUm,2p 0, 


a 


= p G IN 


(78) 


where we do not further need to evaluate V m>a = 2 p here. For the analysis performed in this paper the values 
of the normlization factor (78) for fractional ^ IN is most relevant: 

Explicit determination of (7 7) for the fractional case j ^ IN 


Ifl 


a +1 


sin e “>°+ 


lim 3i(e + i£) 


— Oi—l 


-1 


2 a! sin |a J- 00 


/ OO /TT 

e^|fc|°dfc= - 77^(79) 

-OO < 


a!sin|a v dtf 2 ' 


where sin ^ 7 ^ 0 and we have put a! = T(a + 1). After some manipulations (77) takes the form 


v m , a = (- 1 )" 1 + 1 W .,T.„ . (»a(D - Oa(-I )) 2 " 1 |A|“|a.o, 0 < a < 2m, 

Zj Ct • olxl 2 bt 


a , 

2 ^ 


(80) 
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where this relation is undetermined for integer j 


E 1X1 due to the vanishing of sin = 0 in the nominator 

and the finite difference (Zi>A(l) — ^a(— l)) 2,ri |A| Q! |a=o = 0 in the denominator at these values. However, 
this relation is well determined in the limiting case (in the sense of L’Hopital’s rule) with j = p ± <5 being 
infinitesimally close to integer powers p E I'M o 5 —> 0 (where <5 remains infinitesimally non-zero). This gives the 
possibility to extend (80) analytically also to the “forbidden” integer values. However, we emphasize, that 


there is no need to perform this limiting case if regularization rule (48) is utilized and taking into account 


that the finite difference operator in (80) is cancelling with (52) and hence does not affect the regularized 
representation of the FL (55). 


In (80) the shift operator expression is defined by its binomial series 


(—l) m+1 A 2m (l)|A| a | / \=o = (D x (l) - D x (-l)) 2m |A|“| a =o 

m 9m i 

(D A (i)-D A (-i)) 2m iAru =0 = (-ir E 7 . w ' —ry(—i) p -D A (2p)|A| a | A=0 

p t^ m yrn + p)\(m - p)\ 

2 ml 


(81) 


(ZM1) - D x (-l)) 2m |A| q | a=0 = 2 1 + “(—l) m E 


—J (m + p)\(m — p)\ 


(-I)V 


where the term for p = 0 is vanishing since a > 0 and the symmetry of the terms for has been used. 
The constant (74) yields by using the factorization (75) with (80) and (76) the expression 
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= C 


-i 


7T 2 T(^) 


(-l) m+1 (D X ( 1) - D X (-l)) 2m |An A= o, 0 < a < 2m, ^ W 


m,n,a y -'m,n,a 


2“T(^)T(a + l)sin 


or with (81) 3 


n+1 


Ar, 


27r^r(^) 


E 


(-l) p - 1 (2m)! pa 


m,n,a r (a±n) r ( a + j) sin ^ (jn + p)\{m - p)\ 
For m= 1 (82) takes An=l,n,a = Cm=l,n,a = 2 C n,i- 


(82) 

(83) 


A.2 Explicit determination of the surface integral [/ n>a of (76) 


In this appendix we evaluate the surface integral U n , a dehned by the surface integral (44) on the n-dimensional 
unit-ball. For our convenience we first consider the following integral 


I*,n = [ e - ” 2 |x|“d n r = f |x|“dn(r) ^ [°° e^r^ 1 dr = (E r(^±E (84) 

J J\r\=i Jo 2 Jo 2 2 

where we have put x/r = x. 

Let us assume a € fR and a > — 1: (84) is performed over the entire IR n and we introduce x = rn\ with 
d n r = r ri_1 drdfl(n). Then (84) can also be written as 


I a , n = / e x \x\ a dxH r j =2 / e h dy 3 = 2 tt 2 / e x |x| a dx 

J— 00 J— 00 Jo 

and the latter integral assumes the form of a F-functiori (r = x 2 ) 

[°° -X 2 , ia- 1 f°° -T l/a-l\ 1 .a + 1 

i e Hdl= 2 i e T 2 dT= 2(^-) 1= 2 r( ^- )) 


(85) 


( 86 ) 


We observe that (86) is well defined (non-singular) for a > — 1. and so (84) assumes 
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(87) 



So we obtain with (84) and (87) for the surface integral (76) the explicit form 



( 88 ) 


which is well defined for a > — 1. The dimensions n = 1, 2, 3 are relevant, so we have 



n = 1 


(89) 



n = 3 


n = 2 


(90) 


(91) 
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